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ABSTRACT 

A new technique is described for the full-field determination of soot volume fractions via laser 
extinction measurements. This technique differs from previously reported point-wise methods 
in that a two-dimensional array (Le. image) of data is acquired simultaneously. In this fashion, 
the net data rate is increased, allowing the study of time-dependent phenomena and the 
investigation of spatial and temporal correlations. A telecentric imaging configuration is 
employed to provide depth-invariant magnification and to permit the specification of the 
collection angle for scattered light To improve the threshold measurement sensitivity, a method 
is employed to suppress undesirable coherent imaging effects. A discussion of the tomographic 
inversion process is provided, including the results obtained from numerical simulations. Results 
obtained with this method from an ethylene diffusion flame are shown to be in close agreement 
with those previously obtained by sequential, point-wise interrogation. 


/ 



Introduction 


Of the various optical diagnostics for the characterization of soot volume fractions in flames, methods 
hosed on extinction have been utilized most extensively. 1 ** This is due in part to the relatively large 
optical cross sections, affording the measurement of sman concentrations. The experimental configuration 
Is straightforward to Implement resulting In data that Is simple to analyze. Studies to date have relied on 
sequential point-wise interrogation of the soot field 3 * 4 By translating the position of the optics relative 
to the flame, spatially resolved absorption data is obtained over some region of interest. This represents 
a time consuming process, Inadequate for transient phenomena. 

We have developed a techniq ue for the simultaneous, rapid acquisition of two-dimensional extinction data. 
The method utilizes solid-state detector arrays coupled with collection optics configured to provide depth- 
invariant magnification and characterized angular collection efficiency. Undesirable coherent imaging 
effects and fluctuations in source stability are taken into account to improve measurement sensitivity. A 
system constructed from inexpensive, commercially available components has been used to determine 
fractional absorbances as low as 4 x 10' 3 at a fbll-field acquisition rote of 30 Hz. The system Is simple, 
mechanically stable, and can be configured in a compact fashion, rendering it useful for a variety of 
applications, 

Measurement of Optical Extinction 

Optical extinction can be used to determine soot volume fractions for those particles whose absorption and 
scattering can be approximated by Rayleigh theory. The discussion here essentially follows that of 
Santoro et si in their prior studies of laminar diffusion flames. 3 Neglecting multiple scattering and radiant 
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emission at the detection wavelength, the extinction of optical radiation, by a quasl-continuous medium 
obeys the Bouguer-Lambert-Beer law. 5 

(1) 

ds being the incr ementa l path length and the epical extinction coefficient at the source wavelength. 
Integrating equation (1) along chord-like paths through the flame yields: 

ln(l 1 /i i0 )-/k rt dfl <*> 

In the Rayleigh limit (D < X), the extinction coefficient is given by: 


k Bi ■ (lh) nJ*p (d) D 3 dD 


(3) 


where ffi le the complex refractive Index of the particle, N the number of particles, and E(ffi) is given by: 


BW, ~ Im (^) (4> 

The quantity P(D)dD represents the probability of occurrence of particle diameters between D and D + 

dD, thus Jp ( d ) do taken over all particles is equal to unity. The generalized moment ratio D M is 
defined as: 
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Using this definition, the volume fraction fy can be expressed as: 


(5) 
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From equation* (3) - (6) It follows that: 





<e) 


: ■ 1*1 
Y 67XE(ft) 


( 7 ) 
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For this investigation, a value of ft = A57 ■■ 0,56 ? Is assumed, the exact value Is debated in the 
literature. 2,7 '* Equations <4) and (7) demonstrate that alternative values result In the variation of by a 
multiplicative constant 


Investigations involving both scattering and analysis via transmission electron microscopy suggest that soot 
aggregate structures do not satisfy the criterion of individual spheres with D « and thus the use of 
Rayleigh theory requires validation. 3,9,1 0,11 While tills departure impacts scattering, first-order theory of 
scattering by mass fractal aggregates 12,13 and modeling analyses 14,15,1 * demonstrate that absorption Is only 
weakly affected. Furthermore, the present optical configuration allows the precise specification of 
collection angle, permitting the introduction nf scattering corrections If desired. 


Data Filtering and Tomographic Inversion 

Because optical extinction is path-integrated, a single measurement of transmittance cannot determine local 
values of In addition, the spatial limits over which the absorbing species participate are generally not 

known a prion, precluding the evaluation of the quantity Jds . For axisymmetric flames, the direct 

inversion can be obtai n ed by Abel transformation. 1 7 Numerous approaches to this problem have been 
developed; the reader is lettered to the recent review article by Dasch, 1 * 
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The present work employs the Altered backprojectlon method of RauisrhanUntn and Lakshminaray anan 1 9 
using the algorithm derived by Shepp and Logan. 20 This method has been employed in numerous 
applications involving absorption measurements In flames, due to the resulting accuracy and computational 
efficiency. 3, 21 Before reducing experimental data, the Inversion was validated with numerically 
simulated backprojectlons, including sets which investigated the behavior of the inversion in the presence 
of noise sources specific to coherent lmaglngy\1he two primary noise sources are attributable to the 
coherent effects of speckle and diffraction in the presence of first-order refraction. 

In podm-wlse measurements, the diameter of the detector is normally chosen to exceed that of the 
interrogating beam to reduce sensitivity to beam steering effects. As a result, the measurement is 
integrated over a large number of speckles, reducing intensity variations due to speckle fluctuations. 24 
Any flame perturbs the refractive index field, causing the speckle pattern in the plane of the detector to 
redistribute. The speckle fields are therefore not spatially stationary for reference and data images, but 
appear as additive noise. This Is also true for coherent diffraction effect* from dust particles or localized 
optical defects. Similarly, beam steering causes lateral shifts In the diffraction patterns; the intensity 
perturbations attributable to these shifted patterns are not removed by normalization by the reference 
image. Afi discussed below, Otis effect can be significantly reduced by corrupting the coherence of the 
source. 


Since the derivative-like nature of Abelian inversion accentuates the presence of high frequency noise 

components, 1 ® inversion of data with additive noise is preceded hare by spatial lowpass filtering. Several 

&- 

filtering kernels were investigated, including linearly-weighted convolution and cubic spline smoothing 

i ’* 

(^Simulated data sets were creeled by adding formally distributed random nol6c to numerically generated 
backprojectlons.^ In all cases, the best fit to the original absorbance data was obtained from two- 
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dimensional tensor product B -splines. 25 inversion of tbs filtered numer ics! simulations confirmed the 
accurate retrieval of the original volume fraction distributions. Subsequently, this procedure was 


demonstrated to replicate volume fractions 
laminar ethylene diffusion flames. 3 



Santoro et al for 


Optical Configuration 


The optical configuration is shown in Fig, 1. The beam from a low power laser is followed by a variable 
neutral density filter, then spatially filtered, expanded and collimated. Following the test section are 
Imaging bandpass and neutral density filters; these filters are selected so that no detectable signal results 
from luminous flame emissions. The variable neutral density filter Is adjusted so that the unattenuated 
reference image fills the dynamic range of the detector. The receiving pupil Is located at the front focal 
plane (l.e. Fourier plane) of the weak decollimating lens. This provides a tel acentric configuration 
possessing depth-invariant magnification. 25 In addition, the spat! ally-invariant angular collection efficiency 
Is specified by r^/f, where r p is the radius of the pupil and f is the focal length of the decollimator. The 
virtual image formed by the decollimator is re-imaged onto a monochrome CCD array using a 
conventional oamera lens. The output of the array Is sent to a video-frame digitizer; the results reported 
here were obtained with a bus-mounted eight-bit digitizer residing In a 4B6 personal computer. 

[Hg.l] 

It is essential that the ima ging system be focused on the plane corresponding to the centroid of the first 
moment of the refractive Index taken along the direction of the optical axis. Satisfying this condition 
eliminates intensity variations due to second-order spatial gradients, or shadowgraph terms, 27 which would 
otherwise contaminate the optical extinction measurement Fur radial phase objects, this plane corresponds 
with the axis of symmetry. In practice, optimum focus Is achieved through the use of an axisymmetrlc 
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test object. In this study, a laminar jet of highly refractive gas, freon, was located coaxially with tl» fuel 
nozzle prior to die experiment When die ima ging lens was Improperly focused, a shadowgraph Image 
of the jet wes dearly visible; when correctly focused this image vanished. 


It Is desirable to corrupt the coherence of the laser source to diminish die effects of speckle and 
diffraction. As discussed, these effects are not eliminated by normalization to the reference image, and 
must be attenuated optically. Is the present configuration, this Is accomplished with a rapidly rotating 
glass diffuser plate. Locating this plate near the common focus of the beam expander affords the use of 
a small diameter diffuser, ami minimally docuILUuaies the expanded beam. In using a detector array, the 
divergence of the collimated beam must be equal to or less than dj/(M‘ df), where dp is the lateral 
dimension of a pixel, M is the overall optical magnification, and df is fee diameter of the flame. This 
condition validates inversion by Abel transformation, since non-overlapping chord-like paths are required. 
Figures 2 and 3 depict reference Images obtained with and without the use of the diffuser, fig. 3 clearly 
demonstrates the reduction in speckle and diffraction, in addition, the intensity of die image field is 
spatially more uniform, as the Increased angular spectrum exiting fee diffuser plate effectively o vermis 
the secondary element of the beam expander. 

[Fig. 2,3] 

Results and Discussion 

I m agi n g soot volume fraction measurements were demonstrated using a laminar diffusion flame burner 
similar to that previously reported by Santoro et al. 3 The fuel nozzle has an inside diameter of 10.5 mm 
and is surrounded by a 100 ran diameter annular coflow. The entire assembly is encased in a 405 mm 
long cylindrical brass chimney. Optical access Is provided by two circular ports machined in the wails 
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Of the chimney. These ports are 50 mm In diameter to coincide with the diameter of the expanded laser 
beam. Screens and a ceramic flow restrictor arc placed at die exit of the chimney to achieve a stable 
flame by reducing dr recirculation down the interior aide walls. The fuel port and coflow annulus were 
supplied with ethylene (3.85 cmVsec) and air (713.3 cmVaec), respectively, to coincide with prior 
Investigations. In experiments investigating the effect of temporal averaging, these flow rates were 
reduced to 2 cmVsec and 341.7 cmVsec, respectively, to provide a more fully developed soot shell at the 
same axial locations. 

To account for fluctuations in laser Intensity, each image Is normalized by the average Intensity in a 
specific reference region, here a 10 x 10 pixel area located midway between the flame and the clilmney 
wail. Care must be exercised regarding the conditions under which the reference Image is acquired, in 
this study, the absorption image is captured after the burner reaches steady-state temperature. Best results 
are obtained with a reference image also corresponding to a heated burner, thus the reference Image Is 
acquired Immediately after terminating the fuel flow. The air coflow, however, Is maintained to eliminate 
trBTid^nt refractive index perturbations due to mixing induced by steep thermal gradients at the walls of 
the chimney. These considerations can be appreciated from the observed differences between reference 
corresponding to the rrmritHftnn of a cold or hot burner, or a hot burner with or without the coflow 
maintained . This behavior is attributed to refractive index inbomogendties caused by small geometric 
Imperfections In the burner and buoyant mixing that result In departures from perfect axial symmetry, and 
thus perturb the reference image through second order optical effects. 

[Fig. 4,5] 

Fufl'fleld, single-shot absurbanoe d a t a and the resulting soot volume fraction map for a 3.85 cmVsec 
laminar ethylene diffusion flame are shown in Figs. 4 and 5, respectively. Figures 6 and 7 provide a 
comparison with transmittance data and soot volume fractions at selected axial locations previously 
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reported by Santoro et al. Shown In Hg. 6 is raw tran sm ittance data, as well as the results of filtering by 
linearly weighted, 3x9 pixel convolution and smoothing by tensor B-splines. As is observed in 
numerical simulati ons, the latter provides the hast fit to the original da t a. Slight differences in the 
transmittance data near the burner axis at axial locations of 30 and 40 millimeters above the nozzle tip 
are attri buted to differences in laser wavelength and hence refractive Indices, and small differences in the 
burn er (10.5 millimeter nozzle as opposed to the original 11,1 millimeters, and stainless steel as uppused 
to brass). 

[Fig, 6*7] 

The resulting soot volume fractions shown in Fig. 7 arc in good agreement with previous results, and 
demonstrate shot-to-shot repeatability. The slightly narrower, taller peaks in maximum volume fractions 
are attributed to the two-fold Improvement in spatial resolution (0.1 mm) afforded by the imaging 
configuration. The Increase in noise in the neighborhood of the burner axis Is characteristic of radial 
invention, which to propagate noise from individual backprojections radially Inward. 18 A direct 
comparison Of point- wise and imaging extinction measurements (the latter using our design) performed 
for a laminar ethylene diffusion flame also demonstrates close agreement between soot volume fractions 
obtained by the two methods. 24 

The ability to acquire time-resolved, fuD-fleld data allows the closer Inspection of the effect of time 
averaging, For example, each datum reported in reference (3) represents an average of 100 readings taken 
over a 1 0 second Interval, each reading acquired by lock-in detection operating with a time constant ot 
0.125 seoonds. For the purposes of comparison, to obtain the equivalent of one frame Of 483 x 512 olxel 
image date by point-wise scanning would require an acquisition time of 687 hours, or 28.6 days. In the 
present work, temporal averages are computed fui 2 n image fields, for n ranging In value from 1 to 7. 
For a date acquisition rate of 30 frames per second, a 128 frame average corresponds to an integration 
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time of 4.27 seconds. 


^fwutikcK 

Figure 8 shows ^ercenpbioxbance-d: five axial location* far temporal averages taken over 1, 32, 64, and 
128 fnune averages, Noticeable differences In absorb once levels are observed, while the overall flame 
shape Is preserved. The results suggest that the changes are due to flams movement, because the 
absorbance level Increases at some locations, but decreases at others. Because there is no observable 
direction of convergence a* the averaging time is Increased, we however do not feel that temporal 
averaging will necessarily provide more meaningful data. Rather, we note that, these low-Reynolds- number 
flames are/somewhat-una tabl e, even a ftcr-pfecauttotHsHaksrrto-lsolatc the flameTrom-its-gurToundlogs.' - ' 

r . i ; » i « i 

[Fig. 8] 1 ri )' ^' A;cc p*‘’ r ' $ ; t\ ifkd &><- 

If the difference between the smoothed and original N -point absorbance profiles at a specific axial location 
is used as an indication of the noise level in the original profiles. It Is seen that the 128-frame image- 
averaged profiles are considerably less noisy than the single- shot profiles. The RMS noise, for a 
temporal average of 2n individual measurement Is defined by: 


^RMS = 


\L/2 




( 0 ) 


where A denotes absorbance, and the subscripts o and s denote the original and smoothed profiles, 
respectively. varies from a maximum value of 1.25 (single-shot profile at an axial height of 10 
millimeters), to a minimum value of 0.316 (128 frame-averaged profile at an axial height of 40 
millimeters). Furthermore, it is observed that minimal reduction in N P ^ P occurs beyond a value of n = 
6. For the same axial locations, Fig. 9 shows both uasmoothed and smoothed profiles of single-shot and 
128 frame-averaged data. The agreement between smoothed profiles leads to the conclusion that smoothed 
single-shot data sets provide a reasonable representation of the actual soot volume fractions. 


[Rg.9] 


VI-10 



To account for corrections due to in-sc altering, receiving pupils of different diameters were utilized, 
corresponding to collection angles of 0 . 1 , 0 . 25 , and 0.5 degrees. Under these changes, no appreciable 
differences were observed, a result that Is consistent with Rayleigh theory and the relatively course 8-bit 
digitization scheme employed. Departures from this behavior become difficult to measure in practice, 
since the enhanced forward-scattering resulting from the Increased size of aggregate structures is balanced 
by the larger overall scattering cross sections. The former results In the underestimation of absorption, 
since some photons are forward-scattered into the detector. In contrast, an Increase in scattering cross 
sections accounts for an overall loss of photons from the system, leading to an overestimation of 
absorption. For the 2 cm 3 /sec flames described here, the smallest possible pupil that can be used 
corresponds to a collected angle of 0.086 degrees. For smaller angles, refracted light rays ere intercepted 
by the pupil, resulting In gray-scale schlleren Imaging. The requirement placed on input beam divergence 
to validate Inversion by Abel transformation Implies that for a pupil diameter of less than (dp* f)/(M’ df), 
a collimated source is not required. For the present system, this value corresponds to a collection angle 
of 0.5 degrees. Providing that suitable source radiance Is available to overcome natural flame emissions, 
Imaging absorption experiments could, In principle, be performed using diffuse illumination, e liminating 
coherent effects and providing additional cost savings. This presumes that the source spectrum and 
spectr ally-dependent refractive inde x of soot particles can be suitably characterized. 

Conclusion 

A technique has been developed for rapid acquisition of full-field Optical extinction data. For flames 
possessing suitable geometries, this data can be tomograph! c&lly inverted to provide two dimensional soot 
volume fraction maps. Compared to previously reported sequential, point-wise methods, this technique 
provides a two-fold Improvement in spatial resolution and affords an increase In temporal throughput on 
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the order of 7 x 10 7 . This allows the study of transient phenomena and the investigation of spatial and 
temporal correlations. 

The optical configuration is simple and Inexpensive. While the present system provides the measurement 
of fractional absorbances as low as 4 x 10* 3 tor a spatial resolution of 483 x 512 elements at a temporal 
rate of 30 Image fields per second, further Improvements can be realized through the use of detector arrays 
possessing higher pixel densities and/or more rapid framing rales. Although the sensitivity of the 
instrument Is Improved by corrupting the coherence of the source, these results indicate that the achievable 
threshold sensitivity will ultimately be limited by second-order refractive and coherent optical effects. 

Since coherent optical effects represent an additive noise mechanism not typically encountered in single- 
point measurements, spatial smoothing has been applied tn the resulting full-field data. Numerically 
s i m ulated optical projection data confirms the ability to retrieve absorbance and root volume fraction maps 
In this fashion. The validity of l m a gtn g volume fraction measurements has been confirmed both by direct 
comparison of point-wise and Imaging data, and by the ability to reproduce previously reported soot 
volume fractions in laminar ethylene diffusion fl ames 

The ability to perform temporally resolved measurements has been exploited to Investigate the previously 
unreported effect of temporal averaging. The lack of systematic convergence in the observed soot volume 
fractions leads to the conclusion that, in this application, time-averaging does not necessarily provide more 
meaningful data, but rather that lcw-Reynoids-numher diffusion flames are extremely sensitive to 
perturbation by t h ei r surroundings. Reasonable agreement between spatially-smoothed, single-shot p»^ 
temporally-averaged measurements leads to the further conclusion that the former can be interpreted as 
being representative, on average, of the actual soot volume tractions. 
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Figure Captions 


Fig. 1 Optical configuration for soot volume fraction imaging measurements. 

Fig, 2 Unattenuated reference Image with no diffoser present. 

Fig. 3 Unattenuated reference Image using rotating diffuser plate to corrupt coherence of the Bource. 

Fig. 4 Absorbance map corresponding to the 3.85 cmVsec laminar ethylene UUTuslun flame. 

Fig. 5 Soot volume fraction map corresponding to the 3.85 cmVsec laminar ethylene diffusion flame. 

Fig. 6 Comparison of transmittance data with data of Santoro et al (ref. 3). 

Fig, 7 Comparison of soot volume fractions with data of Santoro et al (ref. 3). 

Fig. 8 Absorbance dam for the 2 cmVsec laminar ethylene diffusion flame demonstrating the effect of 
temporal averaging. 

Fig. 9 Absorbance dam for tbs 2 cmVsec laminar ethylene diffusion flame; unsmootbed and smoothed 
profiles are shown for both single-shot and temporally averaged data sets. 
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Fig. 1 
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Fig. 5 


VI-20 








15 



Fig. 0 
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